Dynamical Scalar Degree of Freedom in Hofava-Lifshitz Gravity 
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We investigate the linear cosmological perturbations of Hofava-Lifshitz gravity in a FRW universe 
without any matter. Our results show that a new gauge invariant dynamical scalar mode emerges, 
due to the gauge transformation under the "foliation-preserving" diffeomorphism and "projectability 
condition", and it can produce a scale invariant power spectrum. In the infrared regime with A = 1, 
the dynamical scalar degree of freedom turns to be a non-dynamical one at the linear order level. 
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Recently, a power-counting renormalizable ultra-violet 
(UV) complete quantum gravity theory was proposed 
by Hofava [H, 0, Sj- This theory is characterized by 
the anisotropic scaling between time and space, so the 
complete diffeomorphism invariance of General Relativ- 
ity (GR) is lost, instead the Hofava-Lifshitz (HL) grav- 
ity is invariant under the so-called "foliation-preserving" 
diffeomorphism. Since this theory was proposed, a great 
deal of efforts have been made, including studies of cos- 
mology lap , I lOL I 111 . Il2l 1 13L I X 51) a nd black hole physics 
[H, lH LL8^ among others THHo, HI HI- Due to the 
differences of diffeomorphism group between HL and GR, 
we expect to see some new dynamical degrees of freedom 
of gravitational fields in HL gravity. Indeed, a new dy- 
namical scalar degree of freedom of gravitational fields is 
firstly discussed by Hofava [l[ in Minkowski spacetime. 
The motivation of this paper is to investigate the dy- 
namical behavior of this scalar mode in a cosmological 
background. 

Let us begin with a brief review about the HL grav- 
ity ■ In order to construct a UV renormalizable quan- 
tum gravity, one possible way is to introduce high or- 
der spatial derivative operators, which make the graviton 
propagator fall off sufficiently rapidly at large momenta. 
On the other hand, in order for the theory to be unitary, 
the Lagrangian can only be quadratic in first time deriva- 
tives of the spatial metric. As a consequence, the theory 
has the anisotropic scaling between space and time. For 
instance, in 3 + 1 dimensions the coordinates (t, x) scale 
as 



t -» t t , 



(1) 



where z is called dynamical critical exponent. In terms 
of ADM formalism the metric can be written as 

ds 2 = -A 2 dt 2 + g VJ (dx i + N^tj (dx j + N j dt^ , (2) 
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where the spatial metric, lapse function and shift vector 
scale as 



9ij -> 9rj 



N -► N , 



A, 



ffZ-l 



A, 



(3) 



The action of the non-relativistic renormalisable grav- 
itational theory proposed by Hofava contains two parts. 
The part of kinetic term is 



Sk = 



dtd 3 x 



IK 



where the extrinsic curvature reads 

Kij = ~ Vi7Vj ~ 



VjA, 



(4) 



(5) 



The part of the potential term in the so-called "detailed- 
balance condition" can be written down as 



Sv — 



dtd 3 x^gN 



k 2 [i e y& 

k 2 h 2 
"8(1 -3/) 



RuV jR l k 



K 2 fi 2 



RijR lJ 



1 - Al 



-R z 



A R - 3A' 



(6) 



where A, k, [i, w, A are coupling constants, e y is the 
antisymmetric tensor defined by £123 = 1 and the Cotton 
tensor reads 



V9 



R j i 



1 



RS j , 



(7) 



The Cotton tensor term in the first line of ©, which 
scales as z = 3, is introduced in order for the theory to 
be power-counting renormalizable in 3 + 1 dimensions. 
The other terms in the second and third lines will make 
the theory undergo a classical flow to z = 1 in the in- 
frared (IR) regime, where the coupling constant A would 
be expected to flow to A = 1. And then GR would be 
expected to be recovered in the IR regime. 

By combining dU) and the full action Si = Sk + S\/ 
can be expressed as 



Si = 



J dtd 3 Xy/g~N ai (K^ K ij -lK 2 ^j 



2 



ijk 



+0xC ij C l > + 11 —R u V J R l k 
V9 

+Ci RijR ij + mR 2 + fi-R + o\ 



with coupling constants 

2 k 2 K 2 fi 

Ql = ^' 01 = ~2^> 71 = 2^ 



(8) 



Ci = — 5— > ??i 



1-4A 
8(1-3A) 4 



o=^4rA, oi = ^^-(-SA 2 ) , (9) 



J1-3A) 



8(1-3A) 



where /i and w 2 are real constants, and have their ori- 
gin as the Newton constant and Chern-Simons coupling 
of Euclideanised three-dimensional topologically massive 
gravity Q. In order to have a real speed of light, A 
must have to be negative for A > 1/3 [3j, which leads 
to a negative effective cosmological constant u\. This is 
not consistent with current cosmological observation. To 
have a positive cosmological constant, one may make an 
analytic continuation of those parameters [f| 



fJL — > tjJL 



w 



(10) 



then the action {8]) changes to 



S 2 = J dtd 3 x«/g~N a 2 (K i:j K tj - IK 2 ) 

+/? 2 CyC ij ' +j2^ F RilV j R l k 
V9 

•C2 /<V>'" + V2R 2 + &R + o 2 



(11) 



with coefficients 

a 2 = ai , (3 2 = -01 , 72 = ~7l j C2 = ~Cl 7 
m = -Vi > 6 = -£1 , cr 2 = -a\ . (12) 

In addition to the gravitational sector {5} or (jlip , we can 
also add some matter sectors to HL theory 



Sm = / d 3 xdty/gN C ma ,tt eI (N,Ni,gij) 



(13) 



Before deriving the constraint and dynamical equa- 
tions of the theory, let us stress that, in order for the 
theory to be tractable, the lapse function N should be 
a function of time coordinate t only, otherwise it would 
lead to difficulties in quantization, at least in the absence 
of extra gauge symmetries. This prescription is known as 
the "projectability condition" [1, 0]- 

Because of the "projectability condition" on the lapse 
function N(t), we can only obtain the spatially integrated 
Hamiltonian constraint by varying the action with re- 
spect to N(t) 

= / d 3 x^{ - a m (K^ - IK 2 ) + f3 m C t3 C l > 



V9 

+ Jn} , 

with m — 1, 2 and 

Jn = -^matter + N- 



RilVjR l k + QrnRijR 1 '' + VmR 2 



SN 



(14) 



(15) 



The equation of motion for Ni gives momentum con- 
straint 



2a m {VjK ji - IV 1 K) + N 







(16) 



The equation of motion for gij is very lengthy, and one 
can find the explicit expression in [8j. 

Now we investigate HL gravity in a flat Fricdmann- 
Robertson- Walker universe 



ds 2 = -dt 2 + a 2 (t)S i jdx i dx j 



(17) 



with the scale factor a(t). For further simplification we 
turn off the matter sector (|13p . In that case, the spatially 
integrated Hamiltonian constraint (TTJ| reads 



H 2 



3a m (l - 3A) 



(18) 



where the Hubble parameter H = a/ a, where a dot de- 
notes the derivative with respect to the cosmic time t, 
while the momentum constraint is trivially satisfied. 
The equation of motion for g^ gives 



2a m (3Z-l) 



H 



-H 



(19) 



By virtue of the above equations, one can easily verify 
that for the action §8§ H 2 < 0, while for the action (TiT|) 
H 2 > 0. Namely in the latter de Sitter solution 

exists. Due to cosmological interest, we will consider the 
latter case and omit the subscript m = 2 in the rest of 
this paper. 

Due to the anisotropic scaling of temporal and spa- 
tial coordinates, the time coordinate t plays a privileged 
role in this theory, and the symmetry of HL theory is 
smaller than the one of GR. Precisely, HL gravity is in- 
variant under so called "foliation-preserving" diffeomor- 
phism, where the coordinate transformations have the 
form 



t -> t = f(t) , x l -> x* = h% x) , 



(20) 



with f(t) is a function of t only. This is a big difference 
from the case of GR. 

In cosmological perturbation theory, the metric pertur- 
bations are usually categorized into three distinct types: 
scalar, vector and tensor perturbations 



Sgoo 



-la 



Sgm — a 2 diB + a 2 



(21) 



a 2 hij - a 2 (d l W J + djWi) - 2a 2 (ipS lJ - didjE) 



3 



where <$>, ip, E, B are four scalar modes, Qi, Wi are two 
vector modes which satisfy d l Qi — d l Wi — 0, and Ziy- 



is the transverse-traceless tensor mode hi 



0. 



Under an infinitesimal "foliation-preserving" coordinate 
transformation 



t^t = t + e°(t) , 
and decomposing e l into 



(22) 



(23) 



with diQ 1 = 0, we can obtain the transformation rules 



<f> -> 4> = <f>- 

B 
E 



# = tp + He , 
B = B-a(Ar), 



E = E 



for scalar modes, 



Wi 



Qi = Qi - a( 



ti 



Wi = Wi + \c. 



2 Si ' 



(24) 
(25) 
(26) 

(27) 



(28) 
(29) 



for vector modes, and the transformation for tensor 
modes is the same as that in GR because of the "foliation- 
preserving" diffcomorphism. Unlike what happens in GR 
Q, however, we can see from (|25l) that we are forbidden 
to choose the "spatially flat gauge" in HL gravity since 
the infinitesimal parameter e° is the function of t only. 
Further, due to the "projectability condition" , the lapse 
function can be set globally to unity, i.e., we can gauge <p> 
mode by choosing a proper initial time. And the residual 
coordinate freedom can be gauged by virtue of (|2T|) . thus 
we can completely fix the coordinates. 

With help of those transformation rules, one can easily 
build up the gauge-invariant variables as 



+ (f>\ 



II EE B - a E 



ip + H I <t>dt' , 

to 



(30) 
(31) 
(32) 



where only two of them are independent, since we have 
$ = {^/H) from (O and (J2U). In the gauge ((f> = 0, 
E = 0), ip coincides with the gauge invariant variable 
"J. Thus it is very convenient to take this gauge in dis- 
cussing cosmological scalar perturbations. In this gauge, 
the metric becomes 

ds 2 = -dt 2 + 2ad l Bdtdx i + a 2 {I - 2tp)5 ij dx i dx j . (33) 



Substituting (|33p into (fTTj) and performing lots of 
straightforward but very lengthy calculations, we obtain 



the action of the scalar perturbations 
S K = 



dtd 3 xaa 3 {(l - 3A) 



2 • 



6Hipip + 3r/r + -ipd B 
a 



-H 2 ip 2 



1 - A 



Bd 4 B 



}■ 



(34) 



3(7 



a 



-a 3 ip 2 



-2£aipd 2 ip} 



(35) 



with d 2 ee S^didj. Although the action (fTTj) contains 
terms such as CijC 1 ^ and RnVjR . the highest order of 
spatial derivatives in the action is d 4 for the sake of the 
antisymmetric tensor e l ^ k and flat universe. However, if 
we abandon the "detailed-balanced condition" , d e terms 
will appear. 

After dropping some surface terms and using the 
background equations (fTS)) (fT9"]) . the spatially integrated 
Hamiltonian constraint becomes 



d 3 x 6a(l - 3X)Hip = . 



(36) 



The equation of motion for diB gives the momentum 
constraint 



<9;{(3A-l)^ + (A-l)-<9 2 s} 



, 



(37) 



and the dynamical equation of motion for the gravita- 
tional field ip is given by 

6a(l - 3A) \$ + "ZHi> +—(d 2 B + 2Hd 2 B 
L 3a V 



12i±M^ = 0. (38) 



This equation is the main result of this paper. 

• When A = 1/3, the first line of (|38[) vanishes, the 
coefficients in second line diverge and the Hamiltonian 
constraint (j3"B"|) is trivially satisfied. So, the evolution 
of ip cannot be determined by the classical equation of 
motion (f3"B"|) . This strange feature indicates that we need 
to take account of the quantum Renormalization Group 
flows at the UV fixed point. 

• When A 7^ 1/3 and 1, substituting the momentum 
constraint (|3"T|) into the dynamical equation (f3"8"|) , we ob- 
tain 



ip + 3Hip~ 



-d ip- 



(1-A)(3C + 8^ V = 



(39) 



a(l-3A)a 2 ' a(l - 3A)a 4 
• When A = 1 , the momentum constraint (|37p becomes 
d4 = , (40) 
with the general solution 

^,x)=/(t)+ 5 (x). (41) 



4 



Furthermore, we can obtain f(t) = const, by plugging 
this solution into the spatially integrated Hamiltonian 
constraint (|36f . As a result, ip is independent of time 
and thus ip is not a dynamical degree of freedom at least 
the linear order level. Of course, it does not exclude ip 
is a dynamical one once higher order perturbations are 
taken into account. 

During the derivation of the above dynamical equa- 
tions, we have not considered the spatially integrated 
Hamiltonian constraint (|36p . Now we explain why we 
can safely ignore this constraint. It turns out to be more 
convenient to discuss this issue in Fourier space 



V>(i,x) 



1 



(2tt) 3 



d 3 k Mty 



then (l36l) becomes 



6a(l - 3A)^ fc=0 = 



(42) 



(43) 



we can see from (|43p that the spatially integrated Hamil- 
tonian constraint only constrains k — mode, while this 
mode can be absorbed into the background. As a result, 
we only need to consider the momentum constraint l|37p 
in HL gravity. 

In what follows we will focus on solving the dynamical 
equation (f3"9")l for the case A ^ 1/3 and 1 . After introduc- 
ing a new field ip = x/a, (|39|) can be written in Fourier 
space as 



Xk( T ) 



a 



\fc 2 B 2 , 2 



(3A-1)A 
(l-A)c 2 



3A- 1 



k 2 - 



■k*T 

2 



X k(r) = 0, (44) 



where prime stands for derivative with respect to con- 
formal time (dr = dt/a) and c for the speed of light 



c = y/ £ / a = -\/K 4 /i 2 A/16(3A — 1). The requirement that 
the speed of light be real implies that A be positive for 
A > 1/3. In the rest of this paper we will concentrate on 
this case. 

Comparing the three terms in the square bracket, we 
can divide our discussions into three cases. 

• Case 1: when k 2 r 2 > T x (2i = A/H 2 ) and k 2 r 2 > 
T 2 (T 2 = \J A/c 2 H 2 ), k 4 term dominates. In this case, 
(1341 reduces to 



Xk( T ) - 



(1 - \) 2 c 2 H 2 
(3A- 1)A 



k 4 T 2 Xk (T) =0 



(45) 



This equation has a general solution with form 



Xfe( T ) = CiParabolicCylinderD 
+ ^ParabolicCylinderD 



-,V2^ 4 kr 



;46) 



parabolic cylinder function. It is easy to check that both 
terms in (05]) are exponentially decaying modes in the re- 
gion —00 < t < 0. In order to see the behavior of this 
solution more clearly, we can use the WKB approxima- 
tion to solve (|45|) in the large k limit. Assuming a trial 
solution in the form of an asymptotic series expansion 



Xk 



(T)=CXpjif>"S„(T)| 
I 6 n=0 J 



(47) 



and plugging it into (|45|) , at the leading order, the solu- 
tion reads 

Xfe(r) = exp |fc 2 C 3 + k 2 ^r 2 } , (48) 

where C 3 is an integration constant and the second term 
in the exponent makes the solution decay exponentially 
when the conformal time r evolves from — 00 to 0, which 
is consistent with the asymptotic behavior of the solution 
(®. 

• Case 2: when k 2 r 2 > T 3 (T 3 = 1/c 2 ) and k 2 r 2 < T u 
k 2 term will dominate over other two terms 



xUr) + c 2 s k 2 X k(r)=0, 



(49) 



where the sound speed c 2 = (1 — A)c 2 / (3A — 1) and a real 
sound speed requires 1/3 < A < 1. The solution of 
is a plane wave solution 



Xk(r) 



cxp{— ikc s r} 
^/2kl^ s 



(50) 



• Case 3: when k 2 r 2 < T 3 and fc 2 r 2 < T 2 , k° term is 
dominant. In this case, one has 



Xk(r) - ^Xk(r) = 



with the solution 



Xk(r) = C A r 2 + ^ , 

T 



(51) 



(52) 



where C4, C5 are two integration constants. We can see 
that, the C4 term is a decaying mode and the C5 term is a 
growing one. For simplicity we neglect the decaying mode 
by simply setting C4 = 0, then we are able to determine 
the absolute value of C5 by matching the absolute value 
of (jSH)) with it when this mode crosses the sound horizon 
(k/aH = -kr = l/c s ) 



|C 5 



1 



(53) 



thus ip will be frozen on the superhorizon scales 



where Ci,C 2 are two integration constants, w = (1 — 
A) 2 c 2 iJ 2 /(3A - 1)A and ParabolicCylinderD[a, x] is the 



(54) 



5 



Finally, by the definition of scalar power spectrum 
■p^(fc) = k 3 \ip\ 2 /2tt 2 , we obtain a scale invariant spec- 
trum 



V* = 



H 2 



\Tt 2 c\ 



(55) 



We can understand these results in the following way. On 
the very small scale (case 1), the fluctuations decay ex- 
ponentially; on the relatively large scale but still deeply 
inside the sound horizon (case 2), the fluctuations oscil- 
late until they crosse the sound horizon; after they crosse 
the sound horizon (case 3), k° term will dominate over 
the other two terms and the fluctuations \ will grow. 
However, this growth is compensated by the growth of 
scale factor, consequently ip is frozen on the superhori- 
zon scale. 

In conclusion, in this paper we investigated the linear 
cosmological perturbations of HL gravity without any 
matter in a flat FRW universe. We studied the gauge 
transformation under the "foliation-preserving" diffco- 



morphism, derived the rigorous equation of motion for 
scalar perturbations and discussed the dynamical behav- 
ior of this mode. Our results showed that this mode 
evolved dynamically with time when 1/3 < A < 1 and 
could produce a scale invariant spectrum. In the regime 
with A = 1, where GR was expected to be recovered, the 
dynamical scalar mode became a non-dynamical one. It 
is of great interest to investigate whether this conclusion 
keeps valid beyond the linear perturbations and some 
matter sectors arc included. 
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